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Abstract: The false vacuum decay in field theory from a coherently oscillating initial 
state is studied for 6 potential. An oscillating bubble solution is obtained. The instanta- 
(]} . neous bubble nucleation rate is calculated. 
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Problems involving quantum mechanical tunneling in a time dependent setting may arise in a 
wide variety of contexts, such as Schwinger vacuum pair production for time-dependent laser 
pulses pQ, pair creation of charged particles in time dependent background electromagnetic 
fields [21 El H], quantum interference in vacuum pair production [5], Hawking radiation from 
black holes [6] , spontaneous nucleation of topological defects in expanding universes [7] and 
false vacuum decay with time dependent initial states or time dependent potentials [U [9] . 

Barrier penetration and tunneling for a particle moving in a one-dimensional potential 
are treated in all textbooks on quantum mechanics. The procedure is by making a WKB 
approximation and expanding the logarithm of the wave function in powers of %. An alter- 
native way to tunneling makes use of the Euclidean-path-integral (EPI) formulation of the 
theory [10]. According to Feynman [11], the amplitude for going from one state to another 
is given by the sum over all paths connecting the states weighted by e lS ^ h , where S is the 
action evaluated along the path. For classically allowed motion, the dominant contribution 
to the path corresponding to the solution of the real-time equation of motion. A convenient 
way to calculate the action is to switch to Euclidean time. In this case, the probability 
amplitude is e~ SE ^ h , where Se is the difference of Euclidean actions between the instanton 
solution (instanton solution: the classical solution of the Euclidean equation of motion with 
appropriate boundary conditions) and false vacuum solution. 
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Most decay of the false-vacuum calculations in single scalar field theory make use of the 
EPI formalism. The Lagrangian of the theory is 

C=±d ll <pd li <j>-V{<i>), 

where V(4>) is a potential which has two nondegenerate minima: <p + (<f>-) is the false (true) 
vacuum. One begins by writing the Euclidean action and the equation of motion. The 
equations are solved to obtain the instanton solution with boundary conditions <p — > ( t ) + f° r 
77—7-00 and (f> ~ for r\ — > 0, where 77 = \/t 2 + r 2 and r is the Euclidean time. The 
instanton solution corresponds to a bubble being nucleated at r = 0. 
The bubble nucleation rate per unit time per unit volume is given by 

T = Ae~ SE/h , 

where Se is the difference of Euclidean action and A is a constant. The relevant solution is 
the one which gives the least action. In flat space and at zero temperature, the dominant 
contribution comes from the unique 0(4)-symmetric solution [T2] . 

As pointed out in [5], EPI has several limitations. We are lost at the outset if <fi couples 
to some external current or field which is time dependent. As an example of this case is a 
scalar field in a Friedmann- Robert son- Walker (FRW) cosmology, since the FRW space is a 
time dependent and cannot be written in static coordinates. Another example arises with 
the theories of two or more coupled field. Also, there is a limitation of EPI formalism in the 
theory of a single scalar field in flat space if the initial field configuration is more complicated 
than simply 4>(x) = <f)+ or time-dependent potential. In this work we can consider the case 
where is homogenous and undergoing coherent oscillations about the false vacuum. 

One approach to overcome these limitations is presented in [Sj. The author studied the 
false-vacuum decay of a scalar field by making use of the functional Schrodinger equation. 
He studied the vacuum decay of a scalar field coupled to a time-dependent external field and 
derived the traversal time for bubble nucleation. 

An alternative approach is presented in [9]. The authors presented a method based on 
WKB approximation combined with complex time path methods, which can be used to 
calculate the relevant tunneling probabilities. They applied their algorithm to production of 
charged particle-antiparticle pairs in a time-dependent electric field and false vacuum decay 
in field theory from a coherently oscillating initial state. For the field theory example, they 
considered the potential discussed in Coleman [TO] . 

K(0) = ^ 2 -a 2 ) 2 + ^(0-a). 

The influence of nontrivial background and decoherence on vacuum tunneling is presented 
in [13]. In this work we follow the algorithm presented in [9], and we discuss the effect 
of coherent oscillating false vacuum state on vacuum decay in the thin-wall approximation 
(TWA), but we choose the 6 potential which was investigated by many authors in the 
context of condensed matter as well as particle physics (see for example [H] [151 US CEO US 
H~9l 120] |2T] [22] 123]). We have noticed that there is a large correction to the nucleation rate 
and the small oscillations about the false vacuum rendered the state more unstable. 

The plan of this paper is as follows. In section 2, the vacuum decay without oscillation 
about the false in TWA is discussed using Coleman's approach. In section 3, decay with 
oscillation about the false vacuum in the TWA is presented based on complex time method. 
In section 4 the structure of the oscillating bubble is obtained, while in section 5 bubble 
nucleation decay rate is calculated. Finally, the results are discussed. 
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Figure 1: The scalar field potential V{4>) with parameters: g = 0.07, A = 2.39, and 5 = 0.2. 

2 Decay without oscillation about the false vacuum: 
Coleman's approach 

Let us consider a scalar field theory with a Lagrangian density 

m = \{W?-v{4>), 

where the potential V((f>) is the effective potential at zero temperature and is given by 

V(<f>) = gcf - 2<A 2 4 + G?A 4 - <5)0 2 , (1) 

we choose g = 0.07 and A = 2.39. The potential is shown in Figure 1, it has two nondegenrate 
minima + (false vacuum) and 0_ (true vacuum) which are all independent of time. 
Let us expand the true vacuum in powers of 5 



To first order in 5, 



(f)_ = A + ei<5 + e 2 <5 2 + ... 



V(cj)_) = -5\ 2 + 0(5 2 )A 



Similarly for the false vacuum 



and 



= + ei 5 + e 2 5 2 + 



for all orders of 5, 
= for all orders of 5 



To calculate the probability of decay of the false vacuum in quantum field theory at zero 
temperature, one should first solve the Euclidean equation of motion of the instanton: 



dV(<f>) 
cUf> ' 



(2) 



with the boundary condition (f) — > <p + as x 2 + r 2 — > oo , where r is the imaginary time. The 
probability of tunnelling per unit time per unit volume is given by 



A e 



S E [ 



(3) 



where Se[<P] is the Euclidean action corresponding to the solution of Eq. fl2]) and given by 
the following expression : 



S f 



d x 



1 /<90\2 

2W 



+ 1(V0) 2 ' + V{4>) 



(4) 



Since we are interested in the lowest-action instanton, we can reduce the problem to one of 
one degree of freedom. If we assume 0(4) rotational symmetry in Euclidean space, then an 
0(4) invariant solution of Eq. ([2]) exists and its action Se[4>] will be lower than that of any 
0(4) noninvariant solution [12]. In this case Eq. ([2]) takes the simpler form 



d 2 <f) 3 dV{<j>) 
drj 2 rj drj d<j) 



(5) 



where rj = y/x 2 + r 2 , with boundary conditions 



<f) — > </) + as 7? — > oo , — = at 77 = 0. 

drj 

We denote the action of this solution by Sq. There is an interesting case (in the sense that 
the action can be calculated analytically) when 



V(<j )+ )-V(^) = X 2 S + O(S 2 )=p 



(6) 



is much smaller than the height of the barrier. This is known as the thin-wall approximation 
(TWA) and the equation of motion (Eq. E]) becomes 



d 2 <f) _ dV(<f>) 
drf d(f> 



(7) 



which can be solved analytically for some potentials. For the 6 potential, the solution has 
the form [TdI EZ] 

A 2 



'wall 



1 + ew 

where /x = y/8g\ 2 , and /z 2 is the second derivative of the potential in the TWA limit evaluated 



at 



The action of the 0(4)-symmetric bubble is equal to 



S 



2ty / df] T] S 

Jo 



--7r 2 p R i + 27T 2 aoR i . 



(9) 
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Here Rq is the radius of the bubble and <To is the bubble wall surface energy (surface tension), 
which is given by 



ctq = / drj 







£d<P ^02(02 _ A 2)2 



v^A 4 
2v/2' 

and the integral should be calculated in the limit po — > 0. 

The bubble radius i? , is calculated by minimizing So, this gives us 

R - 3(T ° 
Po 



(10) 



whence it follows that 



The nucleation rate is then 



T = Ae~ So = ie-V» R o. (12) 



Another parameter which is defined to test the applicability of the TWA is the bubble wall 
thickness L which must be much less than Rq and is given by 

The same results can be obtained using the algorithm proposed in [9]. 
To summarize, in the TWA the instanton takes the following shape: 

{</>+ = A + 0(5), r] « R (True vacuum) 
0waii(r]) = T ^, V~R (14) 
4>- = 0, rj >> R (False vacuum). 

3 Decay with oscillation about the false vacuum: Com- 
plex time method 

In this section we review the results obtained in [9]. We assume that the field is initially 
oscillating around the false vacuum <pf and takes the form 

4>f{t) = 4> + + a sin cut. (15) 

Since the energy is conserved, then E (inside) + E(wal\) of the bubble must equal the 
energy present in the region before nucleation of the bubble: initial- Thus 

47T 

-Ebubbie (inside) = — V(4> t )R 3 , 
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where <p t is the true vacuum and the bubble wall has the energy 



bubble (wall) = 



47r(T bubbl CjR 2 

\ll-R 2 ' 



where 



bubble 



dr 



wall 



^(0bubblc) 2 + ^Obubblc) 2 + ^(0bubblc) 



The initial energy from the false vacuum <f>f(t) has two contributions, namely 



4-7T 

-Einitiai (inside) = — R 3 



47T 



FV D 3 



and 



where 



^initial (wall) 



Ana F /R 2 



dr 



wall 



1 , • 



From conservation of energy 



bubble (inside) + £ bubblc (wall) = initial (inside) + initial (wall) 
which can be written as 

4tt<j e R 2 Att 



1-R 2 



- — Pe R 6 = o, 



with 
and 



„ bubble FV 

<7e=Ve - a E > 



(16) 
(17) 



Pe = p7 -v{<t> t ). 

From the above two equations, we can define the radius of the bubble at some time to as 

3a E 



n 



PE 



and the radius at any later time t (the trajectory) is 



R(t) = V^o + (t - to) 



(18) 



(19) 



The action (S = S^bie — <Sfv) is integrated over an imaginary time contour running 
from some initial time t to t + i7Z , where the bubble shrinks to zero size. The bubble 
action is given by 



bubble 



/ 



dt 



4na L uhhlc (t)R 2 (t)\/l-& + ^-V(<p t )R 3 



(20) 
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where 



bubble 



dr 



wall 



bubble 



1 / • \2 1 / \ 2 

^(^bubblej -^(^bubblcj - ^(^bubble) 



= °E 

while the false vacuum action is 

S F y = -f dt 

where 
and 



wall 



^ r ^bubble 



(21) 



-r = 4 v - 



wall 



dr (ft 2 



FV FV 12 

Pl = Pe -<t>f 



From Eqs. ( J20l) and (f2~Ii) the action is 
S = - fdt 

where 



47T* L (t)R 2 (t)yJl-R*- yPL(t)^ 3 



(22) 



OL(t) =<TE- I dr [Rubble - 0/], 
j wall 

Pl(*) = Pe~ 4>)- 

4 Structure of the Oscillating Bubble 

We calculate the oscillating bubble 0bubbie( r , t) for the 6 potential which interpolates be- 
tween the true vacuum <p t and the false vacuum <pf. Since the initial state oscillates coher- 
ently then it breaks the symmetry of the theory from S0(3, 1) to S0(3). Therefore, Eq. (j2j) 
becomes 



dV 



(23) 



with the potential 

V(<j>) = g<P 6 - 2g\ 2 <p 4 + (gX 4 - 5)<P 2 . 

Following jH], we will find a time-dependent solution 0bubbie( r , t), which will be reduced 
to coherently oscillating field 

<p f (t) = 0+ + a sin cut (24) 

about the false vacuum as r — > oo. The frequency of the oscillations (u) about the false 
vacuum is 



d 2 V 



CO 



< )+ )=2(g\ 4 -6), 



and its range is < u 2 < 4.57. 

We assumed the 0bubbie( r , t) is a function of both space and time and takes the form 



0bubbie(r,t) = <p (r) + a(r) smut. 



(25) 
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After substituting Eq.f l25|) in Eq. (f23|) we get 

fiv 1 

a(r) = 0. (26) 



2 

a"(r) + -a'(r) + 



2 d 2 V 
to - -ttttI 



Now we will solve the above equation of motion in three different regions. 

Firstly, the region outside the bubble (r > R). In this case, O = 0+ (false vacuum), 
0(0o) 0(0+) = w 2 , and Eq. fl26j becomes 

a"(r) + -a'(r) = 0, 
r 

which has a general solution 

C 

a(r) = — + D. 
r 

As r — > oo, ot(r) = 0, hence .D = 0. At r = R, a(R) = C/R = constant which we set it 
equals to a . So, in this region Eq. (J25J) becomes 

0bubbic(^, t) = 0+ + a smut, 

which is the same equation (124)) . 

Secondly, the region inside the bubble (r < R). Again, in this case, 0o = 0- (true 
), 0(0o) ->■ 



vacuum), 77% (0o) — > tOt(0-) = w 2 + A; 2 , where k 2 = 6g A 4 + 185 and Eq. (126]) becomes 



2 

a"(r) + -a'(r) - k 2 a(r) = 0, 
r 

which has a general solution 

, . . sinhfcr 

a(r) = A — . 

kr 

At r = R, a(r) = a and A = a si ^ kR ■ Hence 

R sinh/cr 

a(r) = cto 



r sinhkR 



Note that when a(r) = 0, the oscillation decays to zero inside the bubble. Therefore, the 
thickness of this region (A) is given by 

1 i 

A 



v%A 4 + 185 VTS6 
for small values of 5. Since 

kR~^»l, 

x 2 Vs 

then the solution for a(r) can be approximated to 

R e kr - e~ kr R 



^)=a ^ : kR _:_ kR ^a ^-«y* (27) 



As pointed out in [9], there are three scales which characterizes the structure of the 
oscillating bubbles: the radius of the bubble R ~ 3cr /(A 2 5), the thickness of the bubble 
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wall L~l///~1/ — V(\/8^^ 2 ) an d the thickness of the region inside the bubble 

where the oscillations decay A ~ 1/ (Vl85) and they are related as L « A << R. 

Finally, the region near the wall (r ~ R). In this case 4>o(r) = wa ii(r) and it is computed 
when the potential is degenerate, i.e., when S — > and is satisfying the differential equation 



dr 2 d& 



'wall J 



which has a solution 



0wall(r) 



A 2 



1 + e^ r 

where fi = V"((j)-) = y/8gX 2 up to a correction of first order in 5 and it is the mass of 
excitations around the true vacuum. Since we are working within the frame of the TWA, we 
can neglect the term |a'(r) in Eq. (1261) and we approximate u 2 to u 2 ~ 2gA 4 . Then Eq. (1261) 
becomes 

.2 



a"(r) + 



d 2 V 



'wall 



a(r) = 



which can be simplified to 

a"(r) + 6gX 4 



l + el t(r-R) (l + e n(r-R)y 

By assuming x = fi(r — R), then the above equation becomes 



a{r) = 



a"(x) + 



4 



4 



l + e^ (l + e x 



x\2 



a(x) = 0. 



(28) 



We have solved the above equation numerically which is shown in Figure 2. One can inter- 
polate the solution to an approximate function given by 



a{r) 



B 
7 



0.013O(r - R)) 2 + 5.0 tanh 2 0. 15 (/i(r - R)) - 1.0 



(29) 



Since a(r) — > B in regions r < R, where we know that a(r) = ao, we set B = a^. 

To summarize, we have found a solution for the oscillating bubble in the thin-wall ap- 
proximation 

0bubbie(r,t) = (f> (r) + a(r) smut, (30) 
where (f)o(r) is the static solution given by Eq. ( Tl4l . and a(r) is given by 



a{r) 



a , 



< 20 

4 

a 



0.013(/i(r - R)) 2 + 5.0tanh 2 0.15/i(r - R) - 1 

fl e (r-K)/A 



r > i? + 



<r<i?+f 
r < i?+ f 



(31) 
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Figure 2: The dots represents the numerical solution of Eq. (1281) while the solid line as an 
approximate fit (Eq. (12U|) ). 



5 Bubble Nucleation Decay Rate 

The bubble nucleation rate per unit time per unit volume is given by 

r = Ae -2Im[S(to)l (32) 

We need a time path which shrinks the bubble to zero size. As an example of a path is 



R 2 = nl + (t - t 



which yields to 



t = to + i\JKl - R 2 , for R<TZ C 



We divide the action in Eq. ( 122]) into two parts: 



and 



Si 



to 



to 



dt 



to+iUo 



4na E R 2 ^fl- R 2 -^-p E R 3 



S 2 = - dt 



AnR 2 \ l - R 2 



R- 



rfr (>bubble - </>/ 



2 —R 3 j>t 

3 yt 



(33) 



(34) 



For the first part Si, the calculations proceed as in the static shown in section 2. 

The result is: 



7T 



ImS'i = j^PeK 



(35) 



where p E is the energy density and is given by: 



Pe 



Pe" - Vfa) 
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which is time independent. Note that the oscillation about the false vacuum increases the 
energy density and in the limit a — > 0, pe = Po as expected. 
The surface tension oe is given by 



where a E ubhle is given by 



bubble 



dr 



wall 



„ bubble FV 

a E = <J E - cr E , 



Vbblet/, t)) + - (0 bubblc (r, t)) + V(0 bubb i e ). 



Using 
and 

then 



0bubbie(r,t) = <po{r) + a(r)smut = <p Q (r) + (3(r,t) 
V^bubbic) = V(<M+P(r,t)^(<h) + ^ 2 M)0(0o) 



..bubble 



/ ■ 

J wall 



dr 



± W V(r)cos 2 u;t + l^(r) + \p' 2 (r,t) + <&{r)F(r,t) + Vfa) 



n i x dV , , , 1 n0 . . d 2 V , , . 

= a + O"! + CT 2 , 



where 



and 



cr = / rfr 

■/wall 



^ 2 (r) + \/(0 o ) 



0"! 



wall 



^(r)/3'(r,t) + /3(r,t) — (0 O ) 



which can be shown equals to zero. While 
dr 



a 2 = 



wall 

dr— 

wall 2 

dr— 

wall 2 



i W V(r)sin 2 o;t + ia /2 (r)sin 2 o;t+ l/3 2 (r,t)^(0 o ) + ^ 2 a 2 
^(0 O ) - w 2 ) + « /2 (r)) sin 2 wt + cu 2 a 2 (r) 



a 2 (r 



d<\> 2 



a{r)a"{r) + (a'(r)) 2 j sin 2 wt + u 2 a 2 {r) 
= -tu 2 [ dra\r) = 0.06— La 2 

2 J wall 2 



Rubble = ^ + QM ^L a 2 L 



Therefore, 

while the value surface density due to the false vacuum is 



a E 



dr 



wall 



1 • 



■JM) + U(0/) 



2 2 



i?4 
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Figure 3: The ratio of the relative difference of the radius of the bubble with oscillation 
about the false vacuum and radius without oscillation versus ckq. 



Hence, 



o~e = 0o — 0A7u 2 alL 



(36) 

which is again independent of time, but the oscillation decreases its value and in the 
limit a — > 0, <je = 0o as expected. Moreover, notice that <je equals to zero when 
ctQ = cr /(0A7u 2 L). We would like to see the effect of a on the radius of the bubble. Figure 
3 shows the ratio of the relative difference of the radius of the bubble with oscillation about 
the false vacuum {TV) and radius without oscillation (Rq) versus a^. We notice from the 
figure that at «o = the value of TZ equals to Ro and at «o = 2.57 its value is zero for 
5 = 0.2. So, the allowed value of a is < a < 2.57. 
The second part of the action S2 is given by 



So 



to 



dt 



t +m 



-4:Tioj 2 al 



4irR 2 (t)\/l- R 2 (t) 



wall 



dr{4>l nhhll ... 



(r,t) - $(r,t) --i? 3 (t)^(r,t) 



to 



dt 



to+iRo 



LDR 2 (t)^l - R(t) 2 -\R\t) 



cos cot 



where D = -0.94. Using R(t) = J 11% + (t- t ) 2 and t - t = 1Z zi, then 



So 



-4:7iu! 2 al 



n J° dz (ldkIVi - z 2 - ^nl(i 



— Z 2 



+ iTl f° dzLDKl Vl - z 2 j (e 2i 

+ l n J\z(- 1 -nUi-z 2 )l) 1 -( : 



2io;t0g— 2a;7£oz _|_ g— 2iujto ^IujTZo z 



A7rcu 2 a 2 (h + h + h 



It can be easily shown that 
\mli = 



1 r° 
2 n °Ji dz 



LDTllVT 



\nl{i-z 2 f 
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and 



Jl(0.94)*! + iKj. 



lml 2 = TZ 3 DlJ dzVl - z 2 -(cos2ut ) (e~ 2u)n ° z + e 
= -KlDL^(cos2ut ) J 1 dzVT^e 2 ^ 

2u 4 ZOJlZn v 



-0 

by using Modified Bessel Functions 



u[) y^T(u+l/2)J-i [ t} 6 



and L ~ l/V8gA 2 ~ l/(2w). Similarly for J 3 , 

to) J 1 dz(l-^ 2 ^' 
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1 

Im/ 3 = -^-(cos2u:/i 



— 7r7^o — f cos 2wt ) / 2 (2w7^ )- 



Therefore, 

'0.94 x (2w^o) 2 



ImS 2 = tt^oS (^)(2w72o) + 



32 



/ 94 1 \ \ 

+ {(—)I 1 {2u'R ) + -I 2 {2uj'R )\cos2ut \ (37) 



The total instantaneous bubble nucleation rate is then 



T(to) = Exp 



I 94 1 \ 

+ I (— )/i(2o;^o) + -I 2 (2un ) Jcos2u;to 



(3? 



By fixing the value of 5 to 0.2, we have shown that total action of the instanton (ImS 1 ) 
varies with cto- H has a maximum value (S'max.) at cto ~ 0.12, if we assume cos2c<;to = 1- 
Figure 4 shows a plot of ((S — So)/S max ) versus ao where So is the action give by Eq. ffTTj) . 
At «o — 0, we have S = So while at cto ~ 0.12, most contribution of the action comes from 
the oscillatory part uTZo when it has its maximum value. For a > 0.12, the contribution 
from the oscillatory part starts decreasing and the total action converges to So for higher 
values of «o- So, we conclude that the effect of oscillation about the false vacuum has a 
significant contribution to the tunneling for small specific value of ceo. 

6 Conclusion 

In this paper we have discussed the problem of false vacuum decay in field theory, where the 
initial state consists of coherent field oscillations about the false vacuum for 6 potential. We 
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0.0 0.1 0.2 0.3 0.4 0.5 

Figure 4: The plot of ((S — So)/S max ) versus «o- 

have shown that there is an upper limit for the amplitude of the oscillation of the field about 
the false vacuum. Moreover, The effect of oscillation about the false vacuum has a significant 
contribution to the tunneling for small specific values of the amplitude. The method we have 
used is based on the WKB approximation and the solutions of classical equation of motion 
of the instanton along complex time contour. We obtained a time-dependent decay rate in 
the case of small oscillations. 

The importance of our work is for cosmological models which are based on quantum tun- 
neling, for example: eternal inflation [23], the Hartle-Hakwing-instanton [23], the Hawking- 
Moss instanton [26], the quantum creation of topological defects, e.g. strings and branes 
in a fixed space-time [27]. Moreover, several authors have suggested that string theory in 
four dimensions might have many different vacua [28], which are all represent local minima 
and the tunneling between different local minima is of great importance. Finally, we would 
like to mention here that an important problem which can be investigated is the quantum 
nucleation of cosmic strings and domain walls in an expanding universe. 
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